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MOTION OF THE ANGULAR MOMENTUM VECTOR IN BODY 
COORDINATES FOR TORQUE-FREE DUAL-SPIN SPACECRAFT 


Joseph V. Fedor 


ABSTRACT 


A solution is developed for the motion of the angular momentum vector in body coordinates 
for torque-free, asymmetric dual-spin spacecraft without and, for a special case, with energy dissipa- 
tion on the main spacecraft. Without energy dissipation, two integrals can be obtained from the 
Euler equations of motion. Using the classical method of elimination of variable, the motion 
about the equilibrium points (six for the general case) are derived with these integrals. For small 
nutation an^e, d, the tnijectories about the 9 ■ 0* and d ■ 180* points readily show the require- 
ments for stable motion about these points. Also the conditions needed to eliminate stable motion 
about the 9 ■ 180* point as well as the other undesireable equilibrium points follow directly from 
these equations. These requiremenis are in agreement with C. Hubert (Reference 1). 

For the special case where the angular momentum vector moves about the principal axis which 
contains the momentum wheel, the notion of “free variable’* azimuth angle is used. Physically 
this angle must vary from 0 to 2x in a circular periodic fashion. Expressions are thus obtained for 
the nutation angle in terms of the free variable and other spacecraft parameters. Results for this 
case show that in general there are two separate tr^iectory expressions that govern the motion of 
the angular momentum vector in body coordinates. If the relative angular momentum ratio, h/H, 
satisHes certain inequalities (conditions), one trajectory can be eliminated from consideration. 
Other expressions developed in the paper such as extreme values of the nutation angle in a cycle of 
motion and permissible limits on a dimensionless momentum-energy parameter, H^/2TC, reduces 
to the classical results for spinning torque-free rigid body when the rotor angular momentum is 
set to zero. Energy dissipation on the spacecraft is simulated by a liquid filled ring damper, and a 
simple computer program has been generated that allows the user to determine the motion of the 
angular momentum vector in body coordinates for the special case. Typical trajectories are ex- 
hibited. Also, for small nutation angles, a relationship is expressed between the free variable and 
time, and a damping time constant is derived. 
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MOTION OF THE ANGULAR MOMENTUM VECTOR IN BODY 
COORDINATES FOR TORQUE-FREE DUAL-SPIN SPACECRAFT 

1. INTRODUCTION 

The attitude motion ot dual-spin spacecraft has been examined in considerable detail (Ref- 
erences 1-S, to mention a few) with various mathematical techniques. There appears to be a 
lack, though, of a simple, direct analysis of torque-free, asymmetric duai-spin spacecraft attitude 
motion, unencumbered with complicated mathematical details, that reveals the nature of the 
motion and blends in with the well known classical results of spinning torque-free situation. It 
is the intent of this paper to present such an analysis. Fint, the torque-free situation without 
energy dissipation will be examined and then, for a special case, «nergy dissipation in the form 
of a liquid ring damper on the main spacecraft will be included. 

2. ANALYSIS WITHOUT ENERGY DISSIPATION 

Using a 1-2-3 principal axis coordinate system fixed in the spacecraft and with the origin 
at the center of mass of the spacecraft, the angular momentum of the system, iT, can be written 
as 



where A, B and C are the principal moments of inertia, ccj , coj , ccj are the components of the 
angular velocity vector and h is the relative angular momentum of the wheel which is located on the 

— l|- I -til '4 ^ 'MW 

3 axis. A perfect control system is assumed which keeps the wheel speed constant and hence h. 

The vector equation of motion can be written as 

l^+tlT X TT ■ 0 (1) 

In scalar form. Equation 1 can be written as 

Ad)| + (C — B)cc 2^3 ^ hc»)2 " 0 (2) 

I 


(3) 


B<b2 (A - - hcj| *0 

Clii^ + (B “ A)Ci)|Cj2 * ® 

Multiplying Equttion 2 by o)|, Equation 3 by ^2 and Equation 4 by ci)}, adding and noting that 
the resulting expression is an exact differential which can be integrated gives 

Aw} + Bw| Coi\ • 2T ■ Constant (5) 

Equation S is the rotational energy of the spacecraft with the momentum wheel not rotating and 
it is a constant of the motion. Another constant of the motion, since the motion is torque-free^ 
is the magnitude of the angular momentum vector 

(Awi)* (Bw 2 >^ + (Cwj + h)2 ■ ■ Constant (6) 

Note also, that because of the torque-free nature of the motion, the angular momentum vector is 
fixed in space and can be used as a reference direction to ascertain external viewed motion. 

Referring to Figure 1 which shows a general orientation of a unit angular momentum vector 
in body fixed normalized momentum coordinates 
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one can write 


Acci 

— ■ sin 0 cot ct 


H 

(7) 

or 

■ — sin COS o 

H 

A 

Bu)2 


H 


— ■ sin 0 sin a 

or CC 2 

■ — sin ^ sin a 

(8) 

H 

B 


Cw3 + h 


H cos 9 - h 

(9) 


or «j 


H 

C 



The angle 0 is the conventional nutation angle and a is an azimuth angle. When the normalized 
components of the angular momentum are substit ited in Equation 6, the conservation of angular 
momentum equation, it will be found that this equation is satisfied exactly. This fact is used in 
subsequent development. 

The equilibrium points of the motion will now be obtained using the two integrals of mo- 
tion and the method of elimination of variable. First Acoj will be eliminated which corresponds 
to viewing the motion along the I axis, then Bu >2 which corresponds to viewing the motion along 
the 2 axis and finally t h. Equation 9 is used to write the conservation of angular momen- 
tum as follows 

H2 - (Aw,) 2 + (BW2)2 + H2cos20 (10) 

and the rotational energy equation as 

H2 

2T - Acjf + + (cos d - R,)2 (1 1) 

Where R| « h/H, the relative anguiar momentum ratio of the rotor, cci can be eliminated by 
multiplying Equation 1 1 by A and subtracting it from Equation 10; this gives the following 

H2 - 2AT - (1 - A/BXBw 2>2 + H2 ^cos29 - A (cos 0 - Ri)2j> (12) 

Now the expression in the braces can be written as a perfect square containing cos 0 and some 
residual constant terms. Thus after some rearranging. Equation 1 2 becomes 


3 


0-A/B)(!!^y*(..A/O (co.*-^j.l[R,(..^)4] 03) 


where R 2 U defined as 


J!L 

2TC 


(14) 


As noted previously, Equation 1 3 represents motion projected on a plane normal to the I axis of 
the body fixed coordinate system. It will be noted that the equilibrium point is at 

Ri 


cos $ 


u >2 ” 0 (which implies a ■ 0 or e) and 


R, - — 


I - C/A ' 


(C/A - I) 


(15) 


(16) 


C (C/A - I + R2) 

Now Equation 13 is a displaced conic section locus, and it can be either an ellipse or a hyperbola. 
For an inertia distribution of: A < B < C, the coefficients of the left side squared terms are pos* 
itive and so closed path ellipse trajectories are possible about the equilibrium point for positive 
right hand side, which requires that 


R, > - 


A (C/A - I) 


(17) 


C (C/A - I + Rf) 

In a similar manner for the 2 axis, one can derive the following expression 
0 . b/A) . 0 - B/C, (co, « - ^ [r: (> - -3^) - <>«) 


Again it will be observed, that for this axis the equilibrium point is at 

Ri 


00 s 6 


I -C/B 


(19) 


u)| » 0 (which implies a ■ ±a/2) and 


B (C/B - 1) 
c” (C/B - 1 + R{) 


( 20 ) 
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For th« stated inertia distribution, it can be seen from Equation 1 8 that one coefficient of the 
squared terms on the left side is negative and so the motion about this point is hyperbola in na> 
ture. The slopes of the asymptotes in a plane normal to the 2 axis are given by 

^cos $ - 

These represent separatrices as projected on the 1-3 plane. If the equilibrium point exists, then 
the above implies that no closed path motion will take place about this point. Note further, that 
these separatrices (right hand side) are identical to the simple spinning torque-free case (Reference 6). 

For the two equilibrium points just described, it is interesting to note that compared to the 
sin,tple spinning torque-free case, the presence of a momentum wheel changes the location of 
these equilibrium points only in the Q direction and not in the azimuth (<x) direction. Further, 
as the relative angular momentum of the wheel increases from zero (note Equation 15 and 19) 
the equilibrium points migrate toward the 8 ■ 0” or 180” points, and if the rotor momentum is 
high enough as required by the following 

Ri > I 1 - C/A 1 (22a) 

i> 

and 

R, >11- C/B I (22b) 

the equilibrium points are eliminated (not realizable). This has implications for all attitude angle 
acquisition which will be discussed somewhat more further on in the paper. 

For motion as viewed along the 3 axis. Equation 1 1 is multiplied by C and Equation 10 is 
subtracted from it which eliminates H^cos^d and gives the following intermediate result 

(C/A - 1KAW|)2 + (C/B - 1KBw 2)2 - H2(2RiCOS 8 - R{) - 2TC - H2 (23) 


C/B/ 


± /- 


(B-A) 


A (C - B) 


( 21 ) 
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Now Equation 10 can also be writtrn as 


coftf ■ i (l 


(' 




H2 / 


(24) 


For 9 near zero, cos ^ can be approxinuited by 


cos d W 1 * 

2 


(25a) 


and for 9 near 180*, cos 8 can be approximated by 


. 1 ((Aw,)2+(Bw2)21 

cos 9 W - I + — 

2 H2 


(2Sb) 


Combining Equations 2Sa, b with Equation 23 and rearranging terms, one arrives at 

/AwA^ /BcjA^ 

IC/A .< 1 + R,1 + (C/B - 1 f R| J ) - [l/Rj - (1 - Ri)2l (26) 

for 9 near zero and 

!C/A - 1 - R,1 + (C/B - 1 - RJ . (I + R,)lj (27) 

for 9 near 180*. Equation 26 and 27 represent trajectory motion of the unit angular momentum 
vector if one were to look along the positive or negative 3 axis. The equilibrium points are Ci)| ■ 
0)2 * 0 (which implies 8 ■ 0 or ir for any value of a) and 


1 

R ■ 

(1-Ri)* 


(28) 


for 9 “ 0*, and 


1 

R « 

(1 + R |)2 


(29) 


for 9 ■ 180*. Examining Equations 26 and 27 further, one can conclude that for 9 near zero, 
dosed ellipse like motion is possible if 


C/A - 1+ Ri > 0 


(30a) 
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C/B - I + R, > 0 


(30b) 


and 


R, < 


I 


(I - 


For B near 180*, a similar type of motion is possible if 


(30c) 


and 


C/A - 1 - R, > 0 

(31a) 

C/B - 1 - R, >0 

(3Ib) 

1 



(31c) 

(1 + R|)» ■ 


Note that since T ■ ' , for a given spacecraft moment of inertia and angular momentum, 

the 8 " 0* equilibrium point is a lower rotational energy point than the 8 ■ 180* point. Also, it 
will be noted in Equation 27 that if the roior relative angular momentum Is made laiie enough 
(since physically C/A or C/B is less than 2) so as to make one coefficient negative, no closed path 
motion is possible about the 8 ■ 180* point. 


To emphasize this oonditior, the following is written 


Ri > C/A - 1 

(32a) 

R| > C/B - 1 

(32b) 


which eliminates closed path motion about the 8 ■ 180* point. That is, if inequality 32a or 32b 
is satisHed, then the motion about the 180* point is hyperbolic in nature. If both inequalities 
are satisfied, then closed path motion is still possible provided that 

1 

R > 

(l>Ri)2 

If both inequalities are satisfied in an absolute value sense, the previously mentioned equilibrium 
points about the 1 and 2 axes are eliminated and the only equilibrium points are at 8 ■ 0* and 
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180*. In this case R 2 is bounded by the folio winf 


(1 + R,)» ^ ^ (1 - R,)J 


(33) 


Thus, without energy dissipation, depending upon the magnitude of the rotor relative angular 
momentum, R|, there can be two, four or six equilibrium points for a given dual-spin spacecraft 
moment of inertia configuration. 


It is appropriate at this point to briefly consider energy dissipation. This implies that the 

spacecraft rotational energy (T) is decreasing so that R 2 approaches .. (Equation 33). If 

I U - R|>* 


R 2 were to approach 


, this implies that T is increasing which is contrary to what Ref- 


(I +R|)* 

erence 8 has demonstrated. It is thus seen that the final steady state equilibrium is 8 >0*. 
Summarizing briefly, when the appropriate inequalities are satisfled leaving the 8 ■ G* and the 
180* point u possible equilibrium points, because the 180* point is at a higher rotational energy 
level, this point is also eliminated in a sense that the momentum vector will not move to that 
point with inexorable energy dissipation present. This leaves the 8 ■ 0* as the only remaining 
equilibrium point. The above is important for all attitude acquisition as noted by Reference 1 . 


Another technical tibit can be noted from Equations 26 and 27. If the spacecraft were 
captured about the 8 *• 0* point and it is desired to invert the attitude of the spacecraft, this can 
be done by changing the rotor speed (R}> to a negative value so as to make the 8 ’><' 0* pouit a 
higher energy point and correspondingly make the 8 ■ 180* a lower energy point. The above 
mentioned discussion of energy dissipation would then apply again. 

To add a geometric dimension to the analytic result. Figure 2 shows a sketch of the loci 
of the angular momentum vector on a portion of a momentum ellipsoid for A < B < C with 
three equilibrium points showing (total of six present). 
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Figure 2. Loci of the angular momentum vector for a dual'^pin 
spacecraft on a portion of a momentum ellipsoid for A < B < C 
inertia distribution and h/H > 0. 


3. ANGULAR MOMENTUM VECTOR MOTION ABOUT THE 3 AXIS 

Attention will now be focused on the special but significant case where the angular momen* 
turn vector moves about the principal axis that contains the momentum wheel. This will further 
elucidate and quantify the attitude dynamics of dual-spin spaceaaft. With this motion, the 
azimuth angle that locates the transverse i»mponent of the angular momentum vector must vary 
from 0 to 2ir in a periodic manricr. This angle, a, is called a “free variable'* and is treated as an 
independent variable in what follows. 

As mentioned previously, when the iiormalized components of angular momentum (Equa- 
tions 7, 3 and 9) are substituted into Equation 6, the conservation of angular momentum equa- 
tion, it will be found that the equation is satisfied exactly. Substituting cui, ui* ^3 >nto 
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Equation 5, the enerey equation, and after some manipulation and definition of terms, results in 
a quadratic expression for the nutation angle in terms of a 


G|Gos^ 0 - 2Ri cos ® - Gj “0 


(34) 


where 


or 


and, as before 


Gi 


C C 

Gi ■ 1 - — cos*ot - — sin^a 
^ A B 


I C C 

— - « R? - — cos*o - — sin*o 
Rj ^ A B 


-R| - I 


Rt » h/H, Rz 

living for cos d in Equation 34 gives 


— 

ifc 


cos 6 


R| ± ^ Rj + G| G2 


(35) 


(36) 


(37) 


It wUl be noted that because of the ± sign, there are potentially two trajectories expressions for 
the nutation angle, 6 . The plus sign trajectory is an optional (possible) trajectory depending upon 
R|, the relative angular momentum ratio of the rotor. This is inferred from the inequalities in 


Equations 32a, b. If Equation 32a or b is true, this implies 


the cos 6 equation. The 


x/Rf+C.G, 


> 1 at some azimuth points in 


component will always increase Ri/Gj in the positive 


sign trajectory and thus make the right hand side of Equation 37 even larger than ±1, which, of 

course, is not physically realizeable and thus eliminates that trajectory from consideration. For 

1 

the negative sign trajectory it is of value to let R 2 approach G 2 goes to G^ - 2Ri and 


the cos 6 expression reduces to 


(1 - Ri)2 


COS0 


R, - V (Ri - G,)i 
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which is consistent with R| > C| or G| negative, and gives a value of d ■ 0* as it shou'4. 

The extreme values of the nutation angle in a given cycle of a is of interest and they can be 

de 

obtained by differentiating Equation 34 with respect to o, solving for ^ ; 

do 


dB (C/A - C/B) sin a cos a sin 0 

immm M ^l i i ■ 

da 


Ri - Gi cos 9 


(38) 


de 


and noting what values of a make ~ equal to zero. From Equation 38 it is readily seen that 
this occurs when 


a ■ 0, ir, ±ir/2; for 6 any value. 


dd 

It will also be noted that is zero for 

da 


e • 0, )r; a any value 


which, of course, are the equilibrium poihts. To ascertain if the extreme values are a maximum, 

d^e 

0 ’’^' a minimum, the second derivative must be evaluated at the point. The second derivative, , 
evaluated for sign is 


. d2$ , (C/A -C/B) ^ 

sign . ■ sign — cos 2 a sin 

da2 R|-GiCOS® a-ir/2 


(39) 


As is known, if the sign of the second derivative is negative, then the extreme value in a maximum; 
if the sign is positive, then it is a minimum. Thus, for A < B < C and $ less than 90", the ex* 
treme values, if they exist, are given by 

R, - SORT IRj + (1 - C/BXI/Rj - Rf - C/B)l 


COS 9 


mix 


1 - C/B 


(40) 


and 


cos e 


min 


Ri - SQRT IR; (1 - C/AKI/R 2 « Rf - C/A)l 

TTc/A 


(41) 
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For R| ■ 0, no rotor rotation. Equations 40 and 41 reduce to the ciassicai values of 


BC /I 
1 - C \C ’ 


H2/2T/ 


/ AC /I i~\ 


as it should (Reference 6). It is of interest to note that for cos 0 to be real in Equations 40 and 
41, the quantity in the square root sign must be greater or equal to zero, This requirement can 
be manipulated into the following inequalities: 

^ B (C/B-1) 

^ ^ C (Rj + C/B - I) 


. ^ A (C/A-1) 

^ — 

C (Rj + C/A-1) 

which is consistent with previously derived bounds on R 2 . For A < B < C and for motion about 
the 9 " 0* point, it follows that R 2 is bounded by the following relationship 

B (C/B-1) 1 

C (C/B - 1 + R{) ^ ^ (1 - Ri)2 

Notice that if R| ■ 0, then Equation 44 reduces to the classical result 

B 

-T ^ ^ 1 


HP 

B ^ — ^ C (45) 

2T 

4. ENERGY DISSIPATION ON THE MAIN SPACECRAFT 

The energy dissipation on the main spacecraft caused by a fluid ring damper will now be in- 
cluded in the analysis in a simple manner. Reference 7 developed the optimized energy dissipation 
of the damper per cycle of excitation and it is given by the following expression 
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« 1 ,2 m a* w* 


where 


mass of the fluid 
radius of the ring damper 
amplitude of excitation 


Implied above Is that the change (decrease) of spacecraft rotational energy (Equation 5) is caused 
by the damper. Reference 8 shows this to be valid. If the damper is put on the 2 axis, then 

H 

<u ■ CJ2 * — sin 8 sin ot (A 

* B 

Now a cycle of excitation occurs when a goes from 0 to 2sr. Dividing Equation 46 by 2sr and 
dT 

calling this -r— , we have 
da 


dT -1.2 
da " 17 


m a^ 


Subt^tituting Equation 47 into Equation 48 and symbolically carrying out the integration results in 


T ■ To-0.19 ma7 ~ J sin^dsin^ada 


Where To is the initial rotational energy of the spacecraft. Forming the dimensionless quantity 
Rj - H2/2TC gives 


0.18 R20R3 I sin^® sin 


in^a daj 


Where Rjo is the initial value of R} and 


m a^ C 


a dimensionless quantity containing damper and spacecraft parameters. 
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Now in moit dynamic analyses, the initial iiutation angle is given and not R 20 . Using Equa* 
tion 34, R 30 can be calculated if the initial nutation angle, 6q, initial azimuth angle, Oo, and other 
spacecraft parameters are specified. Thus, 


R 20 


1 

(GjqC0S^9q “ 2R|OOS Oq + R| “ G|o ^ 


(52) 


where 


Gio 



cos^oq - sin^ttg 

B 


In using Equation 52, the maximum value of R^g should be used (which occurs at Oq ■ 0 or ir/2) 
to circumvent square root of nei,ative numbers when calculating the nutation angle. Square root 
of negative numbers indicates that the trajectory is not realizable. 


An interesting aside can be obtained by deriving a small nutation angle damping time con* 
stant for the damper and spacecraft. Using 

cos d » 1 sin^d, 

2 


Equation 34 can be written as 


sin^a 


(I -R|)2 -1/R2 
G, - Ri 


(53) 


Substituting Equation SO for R 2 into Equation 53, differentiating with respect to a to clear the 
integral, canceling where appropriate and averaging over a cycle considering $ quantities constant, 
one obtains 


d« -0.19 R 3 d 

da " (2R, +^/A + C/B - 2) 

Now for $ small, one can physically identify o with time by the following expression 

a »« Xt 


(54) 


( 55 ) 


where t time and X is the body fixed nutation frequency given by 
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(56) 


X ■ SQRT 


' (h-(B-C)c» 3 ) (h -MC - A) 

A B . 


where cjj is considered constant. Equation 56 was obtained from Equations 2 and 3 by consider- 
ing CU| and u )2 varying sinusoidally with Xt. Thus by proper substitutions in Equation 54, a 
damping time constant can be determined 


r ■ 


B2(2R, + C/A + C/B - 2) 


(57) 


0.19 ma2 a 

Note that the numerator of the time constant contains the sum of the left hand side of the in- 
equalities (Equation 30a, b) needed for stable motion about the 6 ■ 0” point. If the numerator 
is close to zero (implying a short time constant) then the expression is not valid because the 
basic assumption of slowly varying quantities is not true. 


5. NUMERICAL RESULTS 

For the special case where the angular momentum vector moves about the principal axis 
that contains the momentum wheel, a simple computer program was generated to calculate nuta- 
tion angle trajectories and other key quantities using the IBM Continuous System Modeling Pro- 
gram (CSMP, Appendix A). The program inputs are: spacecraft moments of inertia. A, B, C, 
relative angular momentum ratio of the rotor, Rj , damping parameter, R 3 , and an initial azimuth 
angle, Oq* and an initial nutation angle, 60 , to specify the initial orientation of the momentum 
vector. The program calculates Rj and picks the trajectory (K « ±1) that will give an initial nu- 
tation angle similar to the input value. The author choose this procedure over the classical 
method specifying R 2 and letting the nutation angle occur where it may. The initial nutation 
angle has more physical meaning to the analyst than R 2 . As previously noted, care must be used 
in choosing Og to circumvent square root of negative numbers. 

The following spacecraft properties were used in the numerical calculations: 

A ■ 91 mass units C/A ■ 1.1538 

B ■ 1(X) mass units C/B ■ 1.05 

C • 105 mass units R 3 « 0.0, 0.1 
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Fifure 3 shows a polar plot of nutation angle versus the free variable a for small motion about the 
$mQ^ point, R| ■ 0.03, without and with damping. Motion is stable, periodic without damping 
u expected, and tends toward zero with energy dissipation. 

Figure 4 shows a Cartesian plot of nutation angle versus the free variable a for R| ■ 0.03 
and R) "0.1 for various initial nutation angles. It will be noticed that some curves go to the 
I) ■ 0" equilibrium point while others go to the d ■ 1 80* point. The local peaks an valleys of 
the curves ore mainly due to the asymmetric moments of inertia. Computer runs showed that 
if 6q were greater than 127*, motion would be toward the 180* point. If do ^*re less than 127*, 
motioi?. was toward the d ■ 0* point. It is interesting to note that d ■ 126.87* is the calculated 
equilibrium point (Equation 19) near the 2 axis. The equilibrium point (d ■ 0* or 180*) is ob> 
viously not unique and the spacecraft can capture right side up or upside down (so to speak) de- 
pending on the initial conditions. 

Figure 5 shows what happens when R| is increased to 0.09. Only the curves for do * 135* 
and 1 70* are shown, since the omitted curves went to the d * 0* point previously. It will be 
observed that both curves now go toward the d ■ 0* point. By increasing R| to 0.09, the in- 
equality of Equation 32b is satisfied and thus, motion about the d ■ 180* point moves hyper- 
bolically away, toward the d * 0* equilibrium point with energy dissipation. Computer simula- 
tions with energy dissipation have shown that when the negative sign trajectory (Equation 37, 

K ■ -1 in the program) is present, invariably the motion is toward the d ■ 0* equilibrium point. 
When the positive sign trajectory is present (K ■ 1) invariably the motion is toward the d « 180* 
equilibrium point. 

6. SUMMARY 

It is felt that the attitude dynamics for torque-free, asymmetric dual-spin spacecraft has 
been developed in a straight forward and elementary manner. Analytic expressions have been 
developed that give results that agree with previously published dual-spin results and that reduce 
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Figure 3. I\)lar plot of nutation angle versus free variable, a. for snuill motion about the point. R| = 0.03, R 3 = 0.0. 0. 1 
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Figure 5. Nutation angle versus free variable, a. for R, = 0.09 and R, - 0. 1 . for 0, - 135*. |70*. 



to classical results for torque-free spinning spacecraft. TImst the neophyte analyst can study 
dual-spin and spinning spacecraft attitude dynamics in a more or less uniform manner. The as* 
sociatcd computer program allows the neophyte analyst to study dual-spin and simple spinning 
attitude motion with and without damping. 
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